Planar spin exchange in LiNi02 
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We study the planar spin exchange couplings in LiNi02 using a perturbative approach. We show 
. . . that the inclusion of the trigonal crystal field splitting at the Oxygen sites leads to the appearance of 

' antiferromagnetic exchange integrals in deviation from the Goodenough-Kanamori-Anderson rules 

' for this 90 degree bond. That gives a microscopic foundation for the recently observed coexistence 

, of ferromagnetic and antiferromagnetic couplings in the orbitally-frustrated state of LiNi02 (F. 

CN . Reynaud et al, Phys. Rev. Lett. 86, 3638). 
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, The compound LiNi02 which was first synthesized in 1958 [Q is known to be a good ionic conductor and therefore 
' suitable as a material for rechargeable batteries. Despite its wide use as electrode material, its electronic and magnetic 
structures are not yet completely understood. Especially intriguing is the absence of any kind of long range magnetic or 
orbital order ^ at low temperature even in the purest samples synthesized up to now. That is especially remarkable 
since the isostructural compound NaNi02 shows orbital order and a collective Jahn- Teller transition at To — 480 K 
from the trigonal high-temperature phase to the monoclinic low-temperature one, followed by an antiferromagnetic 
^ , (AFM) order of ferromagnetic (FM) planes at the Neel temperature of Tn = 20 K |^. To explain the strange 
behavior of LiNi02 the proposal of an orbital liquid was pursued in terms of the SU(4) model 0, §]■ There, 
a symmetry between orbital and spin degrees of freedom is assumed with equal amplitudes for the corresponding 
^ ' coupling terms. In reality, however, the energy scale for orbital interactions is one order of magnitude larger than 
, those for spin exchange interactions as shown by experimental |^ and theoretical studies ||]. This is also indicated 
^ ■ by the difference of To and T/v in NaNi02. More in details, the magnetic susceptibility of LiNi02 shows a transition 
^ , at Tof = 400 K towards an orbitally frustrated state at low temperature. Given an orbital disorder that is frozen in, 
I I ' the magnetic properties at low temperature can be phenomenologically explained assuming FM exchange couplings 
between neighbouring orbitals of different kind {Jdo = (—6.2 ± 2) meV) and AFM couplings between identical ones 
(Jso = (6.9 ± 2) meV) However, that contradicts seemingly the known Goodenough-Kanamori-Anderson (GKA) 
rules [0 which allow only FM couplings for the 90 degree Ni-O-Ni bond in LiNi02 

■ The trigonal (rhombohedral) crystal structure R3m of LiNi02 can be understood by starting from the cubic situation 
T-H , with oxygen and Ni/Li on the sites of the cubes, and with the cubic axes i, y, and z. Perpendicular to the cube 

' diagonal z~x + y + z in fig. 1 one finds alternating planes of Li, Ni, and O. The electronically active Ni02 layer (see 
fig. 1) contains a triangular lattice of the magnetic Ni ions. The ions have a completely filled 2p shell, whereas 
the Ni'^^ ion is in the low-spin electronic configuration (t^g e^) with spin 1/2 

The trigonal distortion changes the bond angle from the ideal value of 90° to 94°, but the six neighbouring oxygens 

■ of Ni stay equivalent. Also the Ni eg orbitals remain degenerate, even in the trigonal crystal field. We now assume 
that we have some kind of orbital disorder that is frozen in for the triangular Ni lattice. The present knowledge 
does not conclusively predict in which way the orbital degeneracy is locally broken. Therefore, we have investigated 

' O ' two situations. The first case deals with orbitals that are oriented along the cubic axes: d^^^p (0) and d^z^^ri (!)• 
The second case chooses two complex combinations of the above orbitals that were predicted in Refs. jll| |lj, |l^ and 
Q ', which have the advantage that they preserve the trigonal symmetry. We will now show that neglecting the crystal field 

■ splitting at the Oxygen sites one finds exclusively FM nearest neighbour exchange couplings. Due to the alternating 
stacking of Li and Ni02 layers, however, one expects a considerable trigonal crystal field splitting of the O 2p orbitals 
(into one doublet and a singlet). That leads to the possibility of AFM exchange integrals in the Ni plane. 
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M ' In the actual calculations we take the Hamiltonian in hole notation. The local part 

H^=H^,+Hl (1) 
contains the Nickel contribution, which for one atom (omitting the site index) is given by 

= ^ ndm + Ud ^ fldm^fldml H ^ ^ ndmndm ^ '^^{2SdmSdm + -jndmfldrn) (2) 

mm m m 

where m and m stand for the two orthogonal Eg orbitals. Ud and Udi are respectively the inter and intra-orbital 
Coulomb repulsions, which are equal when the Hund's intra-atomic exchange Jud is neglected. For an Oxygen atom 
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FIG. 1: Part of the Ni02 layer consisting of the Ni plane (filled circles) and two O planes above (open circles) and below 
(hatched circles) the Ni one. The two coordinate systems used in the calculation are oriented along the cubic {x, y, z) or 
trigonal {x, y, z) axes. Also shown are three orbitals oriented along the cubic 5-axis with the corresponding notation of 
exchange integrals. 



the local Hamiltonian is 

Hp ^ ^ ^n^pn Up ^ ^ ^pnt^pnj, ~t~ ^ ^ ^ ^pn^pn ^ ^ ^ ( ^ Spn Spfi ^" ~'f^p7iflpn) , (3) 

n n n,n n,n 

where n and n are two of the three different 2p orbitals. Up and Upi are the intra and inter orbital Hubbard-like 
repulsions, respectively and J Hp is the Hund's exchange integral. It is well known for p as well as for Cg orbitals, using 
the definitions of Ui, U and Jh as Coulomb's integrals, that one finds the general relation Ui — U ~ 2Jh- However, 
the rotational symmetry at each site for degenerate orbitals demands an Hamiltonian of a more general form than the 
one used in eqs. (§, ||) ||lj, One can keep the truncated parts as in eqs. (||, ^), but then one has to use instead the 
relation Ui = U — Jh to restore, in the case of degenerate orbitals, the invariance of the Hamiltonian under arbitrary 
rotation of the basis at each site. 

The hopping term between Nickel and Oxygen involves the 6 O of the NiOe octahedra. Indexing the Ni site by x, 
the Ni orbital by m, the O site by i, and the O orbital by n, the hybridization part of the Hamiltonian reads 

U^dp = ^ ^ ] ^ ^ (^ai Pina^xma + tix d^xmaPincr) j (4) 



where i™" = (i^")*. The hopping parameters are detailed below. Direct hopping between Ni is neglected. We 
evaluate the spin exchange between two Ni at the fourth-order in perturbation theory in the hopping parameters, 
looking for an effective Hamiltonian of the form 

H'^^^ — Jmm' SxSx' , (5) 

m and m' are the orbital indices of the unpaired holes, x and x' the Nickel site labels. In hole representation, the 
2p levels are empty, while there are three fermions per site in the eg levels. If we neglect Hund's local intra-orbital 
exchange Jnd at Ni, we find in fourth-order perturbation theory the same processes that for the case of one fermion 
at each Ni site. To perform perturbation theory, we make no assumption about the orbital order, looking at all the 
possibilities for m and m', but we suppose that the orbital degree of freedom is frozen in the spin exchange process: 
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TABLE I: Hopping integrals for d^2_y2 and d^i2_^2 
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in the initial and final states the unpaired fermion on Ni is on the same m orbital. Neglecting JhcI-, the result for the 
spin exchange is: 
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(6) 



where the energy needed to transfer a hole from Nickel to Oxygen is £n — — — '2Ud- For a Ni pair, sharing two 
Oxygen, 6 2p levels labelled by n are involved in the intermediate states. The sum on n enables to drop the Oxygen 
site index used in eq. (^), n and n refer to two different orbitals of the same O site. The crystal field splitting at the 
Oxygen sites is taken into account through e„ = A — 5 for two levels per Oxygen, and e„ = A + 5 for the third one. 
The first term of eq. ^) implies a hopping of a hole from the first to the second Nickel; in the second process, the 
two holes jump to the same 2p orbital, or to two different Oxygen sites. The last two terms represent a process where 
the two holes meet on the same Oxygen, but on different orbitals, then the spin flip can act (fourth) or not (third). 
Including Jnd one flnds many additional terms involving spin-flip at the Ni sites. However these terms are typically 
smaller than the previous ones by factors of ^jf^ or ^^fi^. We have explicitely checked that there is no qualitative 
change in the results when we add these terms for Jnd Ud/lO. 

Considering the exchange integrals for the usual Cg orbitals: d^2_y2 (0) and d^^2_j,2 (1) it can be seen from fig. 1 
that not only the two Ni orbitals need to be specified, but also the particular Ni pair: this leads to 6 different spin 
exchange terms J^j^, Jqq, J(li, Jqi, Jj^, jj'^^. The two degenerate 2p orbitals are labelled px and py, the last one pz- In 
terms of the cubic Px, Py, and p^ orbitals they read 



Px = -^{px-Py), 
Py = ~ 
Pz = -^(Pi +Py+Pi) 



(7a) 
(7b) 
(7c) 



The hopping parameters between the Nickel orbitals (0) and (1) and these Oxygen orbitals, depending on the O atom 
involved, can be obtained by linear combinations of the hopping integrals between the Ni orbitals and (j)x,PyTPz) 
which are specified in table I. An arabic index enables to distinguish between all the 6 Oxygen of the NiOg octahedra, 
in correspondence with fig. 1. 

The parameters Ud = 9.5 eV, Up = 4.6 eV, and io = 1-3 eV were taken from Ref. The charge transfer energy 
A = 4 eV was roughly estimated from the difference of the centers of gravity of the 3c? and 2p levels in a bandstructure 
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FIG. 2: Exchange couplings between orbitals oriented along one of the cubic axes 



calculation which allows also to estimate the crystal field splitting S from the bandwidth of the upper 3<i Cg band. 
Namely, for 6 — and only nearest neighbour hopping in the tight-binding Hamiltonian (eq. (^), the bandwidth 
would be zero. We obtained a value for S of about 1 cV, but we keep S here as a parameter. The result for the 
different spin exchange couplings as a function of the crystal field splitting S is presented in fig. 2. For 6 = 0, the only 

contribution to exchange is from Jnp {Jhp = Up/IQ). We used Upi = Up — Jhp (see above). A splitting of |^| « 0.5 

eV can change the ferromagnetic character into a antiferromagnetic one for the exchange integrals Jqq, Jf\, and 
(The exchange coupling Jqq is identical to zero in the present approximation.) The exchange terms between different 

orbitals J^^ and Jq^ stay ferromagnetic (see Ref. [||). 

It seems as if an antiferromagnetic exchange integral between equal orbitals is in contradiction with the existence 
of ferromagnetic planes in NaNi02. One should keep in mind however that the low temperature phase of NaNi02 is 
not trigonal but has a large monoclinic distortion which leads to a different crystal field splitting at Oxygen. One 
could speculate that the monoclinic distortion in this coumpound can be described in our calculations by a reduced 
effective parameter 5ef / shifting all the exchange integrals into the ferrromagnetic region. 

At present there is no experimental indication for macroscopic Jahn- Teller distortions of the trigonal crystal struc- 
ture of LiNi02 even at very low temperature. Therefore it is tempting to search for an orbitally-frustrated state in 
terms of orbitals that obey the trigonal crystal symmetry. For that purpose, noting that the Cg stay degenerate, we 
can use instead of dx2-y2 and d3j2_^2 any linear and orthogonal combination of these two. The particular choice 

da = —^{di:2_y2+id:i2^_r^), (8a) 

d/S = --i=(d£2_j;2 - id3j2_r2) (8b) 

corresponds to orbitals which satisfy the trigonal symmetry: under a rotation of ^ around the axis perpendicular 
to the Ni plane (fig. 1), these states are just changed by a phase factor. These complex orbitals were previoulsy 
proposed for manganites 0, |l3|. It cannot be a priori excluded that also in LiNi02 one of these complex orbitals 
is preferred at each Ni site. 

In terms of these new orbitals the calculations get simpler: there are only two different exchange integrals: Jg = 
Jaa = <//3/3 and Jd = Jaf3, and there is no need to consider different Ni pairs, due to the trigonal symmetry. However 
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FIG. 3: Exchange couplings between orbitals oriented along the trigonal axis. 

it is not possible to establish a simple link between these two new exchange terms, and the 6 previously studied: this 
is due to freezing the orbital degree of freedom in the calculation of the spin exchange. Using the transformation (H) 
and the table (Q) one can easily specify the hopping for complex orbitals to be inserted in eq. (||). 

The results are presented in fig. 3. At the same critical S- value as in fig. 2, the sign of Jd changes to an AFM 
coupling for large 6. But now, the effect is reversed in comparison to fig. 2 or Ref. it is the exchange between 
different orbitals that becomes antiferromagnetic. 

In the basis of complex orbitals, we can also consider the pseudospin exchange. Due to pseudospin anisotropy, there 
are many terms in the effective Hamiltonian involving the pseudospin operators. These terms are usually of the order 
of Up4/{A^{2A + Up)) for 5 = 0. This confirms that they are larger than the spin exchange terms characterized by 
J Hp- We have considered more in details the z component of the pseudospin exchange which corresponds to a term 
j'^TfTj for a pair of Ni sites labelled by i and j, with = +(— )| depending of the orbital involved da {dp). We 
suppose that the spin degree of freedom is frozen at each site. Somewhat surprisingly when the Hund's exchange term 
for Oxygen, as well as as the splitting of the 2p levels are neglected, we obtain that Jj — 0, indepently of the spins. 
That shows the highly anisotropic character of the pseudospin exchange. But one should note that this result is not in 
contradiction with the results from Mostovoy and Khomskii (who found a positive value for the orbital exchange) 
since they analysed the usual d^2_y2 and d^^^-^^ orbitals which are eigenstates of the pseudospin operators in the 
present basis. The pseudospin-rotational symmetry being not required as for spin, the result jJ = does not rule out 
a pseudospin exchange along the others directions. The evaluation of other pseudospin-pseudospin or pseudospin-spin 
coupling terms will be detailed in a forthcoming paper using an alternative method with a direct diagonalization of 
the hopping part of the Hamiltonian which enables to evaluate the different terms in a more direct way. 

In conclusion, we have shown that a proper inclusion of the trigonal crystal field splitting at the Oxygen sites of 
LiNi02 leads to a coexistence of AFM and FM couplings. Orientating the orbitals along the cubic axes we find AFM 
couplings between equal orbitals and FM couplings between different ones, which is at least in qualitative agreement 
with the analysis of Ref. |^. Investigating the orbital symmetry breaking in LiNi02 in terms of complex orbitals is 
fascinating from the point of view of symmetry: one has only 2 orbitals and 2 different exchange couplings. Also in 
that case we found a coexistence of FM and AFM couplings which may explain the absence of magnetic long range 
order at low temperature in LiNi02. But the actual signs of exchange between different and equal orbitals are just 
reversed in comparison with the analysis of Ref. . 
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